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Abstract 

Generalising expressions given by Komar, we give precise definitions of gravi- 
tational mass and solitonic NUT charge and we apply these to the description of 
a class of Minkowski-signature multi-Taub-NUT solutions without rod singulari- 
ties. A Wick rotation then yields the corresponding class of Euclidean-signature 
gravitational multi-instantons. 
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1 Introduction 



In many respects, the Taub-NUT solution [1] appears to be dual to the Schwarzschild 
solution in a fashion similar to the way a magnetic monopole is the dual of an electric 
charge in Maxwell theory. The Taub-NUT space-time admits closed time-like geodesies 
[2] and, moreover, its analytic extension beyond the horizon turns out to be non Hausdorff 
[3]. The horizon covers an orbifold singularity which is homeomorphic to a two-sphere, 
although the Riemann tensor is bounded in its vicinity. These pathologies lead to the 
view that the Taub-NUT solution is not physical. 

Nonetheless, despite the fact that no magnetic monopole has yet been found in our 
universe, such magnetic dual solutions play an important role in quantum electrodynam- 
ics and especially in its non-abelian generalisation, namely Yang-Mills theory. Moreover, 
the stationary solutions of the Maxwell-Einstein equations admit a non-linearly realised 
SU (2,1) symmetry group [4] which generalises the Ehlers group and which mixes to- 
gether the electromagnetic and the gravity degrees of freedom. This generalises to a 
large class of theories, and in particular to ones that can be embedded into supergravity 
theories. Despite the fact that this has not been proven so far, these symmetry groups are 
believed to act on the non- stationary solutions as well. The major difficulty in formulat- 
ing such symmetries comes from the fact that Einstein's theory is highly non-linear and 
consequently its dualities are poorly understood beyond the linearised level. The aim of 
this letter is to understand more closely the duality relations within Einstein theory by 
exhibiting their similarities with the example of Maxwell theory, and more specifically 
the similarities between NUT sources and magnetic monopoles. 

While a magnetic charge can be expressed in terms of a current associated to a vector 
field dual to the standard Maxwell potential, its expression as a function of the standard 
vector potential corresponds to a topological invariant of the associated fibre-bundle 
geometry. In this letter, we define the NUT charge, in a similar way, as a topological 
invariant associated to time-like three-cycles. We also generalise the Komar mass to the 
case where there is no space-like slice with compact boundary in the asymptotic region. 
These definitions involve a fibre-bundle construction which is very reminiscent of the 
one appearing in Maxwell theory. In this case the U{1) fibres are orbits of the time-like 
isometry. 

We exhibit the similarities between the Komar NUT charge and magnetic charge 
through a consideration of explicit solutions involving several NUT sources. Indeed, we 
will give an infinite set of new regular solutions of the Einstein equations with an arbitrary 
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odd number of NUT sources. We obtain these by acting with the U{1) duahty group for 
stationary solutions on muhi-black hole solutions with both negative and positive masses. 
We define a coordinate patch that permits us to avoid Dirac-Misner string singularities. 
Then we show how one can avoid the conical singularities usually appearing in multi- 
black hole solutions by choosing adjacent NUT charges to be opposite in sign. Since both 
negative and positive NUT charge singularities are covered by horizons, these solutions 
define space-times which have no more (albeit also no less) pathology than the ordinary 
Taub-NUT space-time. 

The resolution of the Dirac-Misner string singularities requires the quantisation of 
NUT charge. The Chern class of an associated fibre-bundle geometry is understood 
to count the relative number of fundamental NUT charges of a given spacetime. The 
timelike three-cycles surrounding several NUT charges turn out to be diffeomorphic to 
Lens spaces L{\N\, 1) = S'^/Z|jv|, where N is the relative number of fundamental NUT 
charges that lie inside the interior of the corresponding three-cycle. 

We discuss in a final section the Euclidean analogues of these multi-NUT solutions 
which are slight generalisations of the instantons described in [5] 

2 Komar NUT charge 

A. Komar defined the mass for asymptotically Minkowski solutions of the Einstein equa- 
tions through an integral over the boundary of an asymptotically space-like hypersurface 
V in spacetime [6]. Given an asymptotically Killing time-like vector k, = K^d^, the metric 
permits one to define the 1-form g{K) = g^^n^dx'^ , and the Komar mass is then given as 
a function of the 2- form K = dg{K) by^ 



where -k is the Hodge star operator. Comparing then this formula to the ones defining 
the electric and the magnetic charges 



"'^Komar proved in [7] that if k is chosen to be orthogonal to a family of minimal hypcrsurfaces, 
then the Komar mass will be positive if V is chosen to be one of these hypcrsurfaces. However, Misner 
then showed in [8] that this prescription is either inconsistent or impossible to achieve in some relevant 
examples. Here we will not insist on this orthogonality prescription and the Komar mass consequently 
will not be necessarily positive. 




(1) 
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it seems natural to define the dual mass as the integral 



n = — 




K . 



(3) 



However this integral is trivially zero because of Stokes theorem, as also would naively 
be the one defining magnetic charge. Nevertheless, not all asymptotically flat space- 
times admit a well-defined asymptotically space-like hypersurface. This is the case for 
instance for the Taub-NUT space-time, for which the r = const slices of any space-like 
hypersurfaces are not closed manifolds [2]. 

Let M be an asymptotically fiat space-time. Strictly speaking, we assume that M 
admits a function r which goes to infinity at spatial infinity and which defines a proper 
distance in this limit, g^'^d^rd^r 1, and we assume that all the components of the 
Riemann tensor in any vierbein frame go to zero as 0{r~^) as r ^ +oo. We consider 
stationary solutions; k is then a Killing vector and the second-order components of the 
Einstein equations can be written as 



which is very similar to the Maxwell equation. We choose the function r to be invariant 
under the action of the time-like isometry, and choose the squared norm of the time-like 
Killing vector gp^pK^K!^ = —H to tend to —1 -|- 0{r~^) as r — -|-oo. We assume that 
the action of the time-like isometry is free and proper on the domain of M where the 
function H is positively defined. This implies that M admits an Abelian principal bundle 
structure over a Riemannian three-fold V on this domain. 

If this principal bundle is trivial, it admits a global section s which defines an em- 
bedding of V into a space- like hypersurface of M (outside the zeros of H) . Otherwise it 
only admits a patch of local sections defined on each open set of an atlas of V, which 
we denote collectively by s as well. A natural generalisation of the Komar mass formula 
thus consists in defining it as the integral of the pull back s* of the 2-form -kK 
over dV. In order for this integral not to depend on the local trivialisation, -kK must 
be horizontal and invariant in the asymptotic region. It is trivially invariant since it is 
built from the metric and the Killing vector, and the horizontality condition is given by 
asymptotic hypersurface orthogonality, i.e. ii^-k K ^ as r ^ +00. Because of equation 
(4), dif^-k K = in the vacuum and if space-time is simply-connected, there exists a 
function B such that ii^-kK = dB. The horizontality condition for -kK in the asjTiiptotic 
region is then equivalent to the fact that B tends to zero as 0{r~^) as r ^ +00. The 



dkK = 2kdx''RppK'' = 16nGkdx''{Tpp - \g^,uT)K- 



V 



(4) 
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2-form K is also trivially invariant and its horizontality condition ^^^K = dH — is 
satisfied because the function H tends to unity as r — > +00. 

We accordingly define the mass m and its dual, the NUT charge n, by the following 
integrals 

m = — [ s*^K n = —[s*K. (5) 



OTT Jqv stt Jqy 

By construction, the 1-form g{K) is invariant under the action of the time-like isom- 
etry, and since ingii^) = —H — > — 1 as r — ^ +00, it defines a connection on the principal 
bundle in the asymptotic region. The NUT charge is proportional to the Chern class 
of the principal bundle over dV , and is thus non-zero only in the case where the latter 
is non-trivial. Real line bundles over a compact surface always have a vanishing Chern 
class, and a non-zero NUT charge implies therefore that time-like orbits are compact. 

One defines electric and magnetic charges in the same way by requiring both the 
Maxwell potential and its dual to be invariant under the covariant action of the time-like 
isometry in the asymptotic region, i.e. ii^F i^^-k F (9(r~^). The whole construction 
can be generalised to non-stationary space-times, as long as C^g^y tends sufficiently fast 
to zero as r — s> -|-oo. 

Let us now express the mass and its dual in a more explicit way. We choose coordinates 
for which = in such a way that the metric is given as follows 

ds^ = -H{dt + Bidx'Y + H~^%j dx'dx^ . (6) 

The vacuum Einstein equations then give di^,-^ K = and 

i^i^K = -H^^Ei^^djBk dx' = dx'd.B (7) 

where Latin indices are raised and lowered with the three-dimensional metric 'jij. The 
asymptotic horizontality conditions for K and -kK are satisfied if 

H = l-'^ + 0{r-') B = -y + 0{r-') . (8) 
Then icK and K have the following behaviour in the asymptotic region 

~ ^Eij^H-^dkHdx' A dx^ K ~ -diI3j dx' A dx^ . (9) 

If we assume furthermore that is asymptotically Euclidean, one may verify that the 
parameters m and n appearing in (8) are truly the mass and NUT charges defined by 
(5). 
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One obtains s*K = 2nsmdd9 A dip in polar coordinates on dV = S"^. Bi is thus only 
globally defined up to a constant, and the time coordinate also is not globally defined over 
the two-sphere. One defines t+ and t_ on the north and the south pole of the two-sphere 
respectively. These coordinates are related by 

t+ = t_ - Arvp . (10) 

Since </:> is a periodic coordinate of period 27r, the time coordinates t± must be periodic 
of period Svmo, such that = ^ is an integer. The integer N parametrizes the Chern 
class of the principal bundle over S*^, and the r = const slices for r sufficiently large are 
diffeomorphic to the Lens space S'^ /'L\n\- 

The Lens spaces are usually studied as Riemannian three-folds, but they also admit a 
pseudo-Riemannian metric, as does any U{1) principal bundle over a Riemannian mani- 
fold. If we define a connection uo on to the principal bundle, as well as the pull-back of 
the metric 7 on the base by the bundle projection vr, then uj®uj + tt*'-) gives a natural Rie- 
mannian metric on the principal bundle, and — u;®u; + 7r*7 a natural pseudo-Riemannian 
metric. 

S. Ramaswamy and A. Sen obtained a similar result in [9], where they defined the 
NUT charge as a dual of the Bondi mass instead of the Komar mass. The Bondi mass 
and its dual are defined using integrals involving respectively the Weyl tensor and its 
Hodge dual. 

The f/(l) principal bundles over 5*^ are classified by their first Chern class, which is 
unity in the case of the Hopf fibration of S'^. By analogy with the case of the Maxwell 
theory for which the Chern class determines the relative number of fundamental Dirac 
monopoles, we will wish to interpret this integer as the relative number of fundamental 
NUT sources in General Relativity. This interpretation turns out to be right, as we shall 
see in the following. 

3 Multi-Taub-NUT solutions 

We now want to consider axisymmetric stationary solutions of the Einstein equations 
with several NUT sources on the axial symmetry axis. We use Weyl coordinates in which 

ds^ = H-^e^"{dz^ + dp^) + p^H'^d^"^ - H{dt + I3d^f . (11) 
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For example, in the case of the Taub-NUT solution of mass m and NUT charge n, the 
Weyl coordinates are related to the Schwarzschild ones by^ 

p = ^r"^ — 2mr — sin 6 z = {r — 'm)cos6, (12) 

in terms of which the metric is 

ds'^ = -H{dt± + 2n{±l - cos d)dipy + H-^dr'^ + (f^ + n^) [dO^ + sin^ edif"^) (13) 

with 

H='l^pL^. (14) 

In Weyl coordinates, equation (7) reduces to 

p^^H'^dpB = -d,B p-^H'^d,B = dpB (15) 

and B is the imaginary part of the so-called Ernst potential, £ = H + iB. This latter 
satisfies the Ernst equation 



{£ + £*) (^9/ + dp^ + ^dp^ £ = 2d,£d,£ + 2dp£dp£ . 



(16) 



For static solutions, the Ernst potential is real and the Ernst equation reduces to the 
linear differential equation 

d,^ + dp^ + ^dp^\n£ = . (17) 

The product of several real Ernst potentials thus gives a new solution. This permits one 
to obtain the Ernst potential of multi-black holes solutions as 

where 2ri = + rt^ with 

r,± = ^{z -Zi± q)2 + p2 (19) 

and where Zi and q define respectively the position and the (possibly negative) mass of 
each of the h black holes. When all masses are positive, these solutions are always known 



^Note that the radius r that is commonly introduced in Weyl coordinates is not the Schwarzschild 
radius f, but is related to it by r = f — to. 
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to suffer from conical singularities unless one considers an infinite chain of black holes 
[10]. 

A nice way to interpret the NUT charge as a dual mass comes from the fact that 
the stationary solutions of Einstein's equations admit a nonlinearly realised U{1) Ehlers 
symmetry [11] which rotates the mass into the NUT charge in the case of the Taub-NUT 
solutions. This U{1) acts trivially on the conformal factor a and modifies the Ernst 
potential as follows 

^ , , cos aS — isina , . 

S(a) = . 20 

cos a — z sm a c 

Acting this way on the Ernst potential (18), one gets 

^ _ cos a]\{ri - Cj) -ism.aY[{ri + q) ^^^^ 
cos a]\{ri + d) - isma]\{ri - q) 

where e^*"cj = mj + irii. We then derive the potentials for the metric^ 

H= ^ ^^[f 5 = 6-2Vn,^ (22) 

cos'^ a [[[ri + CiY + sm. a [[[Ti - CiY 

where h is an undetermined integration constant coming from the duality relation (15). 
Note that the potential is a sum of potentials for ordinary Taub-NUT solutions indi- 
vidually centred at Zi. There is one horizon on each segment p = 0, — |cj| < z < Zi + \ci\. 
Let us consider that they are all separated, i.e. that 

Zi-i + |ci_i| < Zi - \ci\ . (23) 

Between each adjacent pair of horizons, there is a Dirac-Misner string singularity related 
to the fact that the 1-form dip diverges on the symmetry axis p = 0. The Dirac-Misner 
string singularities are located on h + \ segments Dj, on which p = and Zi-i + |ci_i| < 
z < Zi — \ci\, where we understand — oo < z < Zi — \ci\ and Zh + \ch\ < z < +00 for Di 
and -D/i+i respectively. In order to avoid such a singularity, the potential 13 must vanish 
identically on each of these segments. On the segment Di, rj = z — Zj for j < i and 
Vj = —z + Zj for j > i, so one has 

i-l h 

=6,-2 J]n,+2^n, =0 . (24) 
j=i j=i 



To derive the formula for B wc observe that H ^dB ~ — 2 V "Z'^''',) , and wc make use of the 
identities + = 1 and ± = ± q to show that 

pdyJl _ „ Z- Zj pdpTj _ Z- Zj 

2 2 P 1 2 — 

C-i 7*2 Ti Ci T*j 
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Exactly in the same way as for the ordinary Taub-NUT solution [2], in order to avoid 
Dirac-Misner string singularities, one must define h + 1 open sets f/j, such that u'll^^Ui 
covers space-time outside the horizons. We define each Ui as the complement of the 
domain Uj^iDj in M. At the intersection between Ui and t/j+i, the corresponding time 
coordinates are related by 

ti+i = ti- Auiip (25) 

and B is given by 



-I h h 

z — z~ 



5|^. = 2^n,-25^n,-2^n,— (26) 

j=i j=i j=i 3 

on Ui, in such a way that dti + B\u-dip is globally defined on M. 

Since if is a periodic coordinate, ip ~ ip + 27r, consistency requires the time coordinate 
also to be periodic, that is tj ^ tj + Snni for all rii. In order for the manifold to be well 
defined, all the NUT charges Ui must thus be integral multiples of a given fundamental 
charge no, so tj ^ tj + STrno- 

We thus conclude that, even on a purely classical level, the existence of more than 
one NUT charge on a manifold implies the quantisation of these charges. In fact, this 
quantisation already occurs in Maxwell theory if one considers that its solutions are the 
connections of U{1) principal bundles over space-time for which the curvature verifies 
the equation d-k F = 0. Indeed, the global definition of the Maxwell connection on the 
principal bundle similarly requires all magnetic charges to be integral multiples of a given 
fundamental charge. 

The one-form cu = ■^{dti + B\u.d(p) defines a connection on the U{1) principal bundle 
over V. For any two-cycle of V surrounding a subset / of the NUT charges, one computes 
that associated Chern class to be 



iV, = -^n, . (27) 



The time-like three-folds that surround the NUT charges within / are thus diffeomorphic 
to the quotient of by Z|Ar^| acting as a discrete subgroup of U{1), yielding a Lens 
space. We thus interpret the Chern class Nj of a two-cycle as the relative number of 
fundamental NUT charges inside its interior. 

As for multi-black hole solutions, the multi-NUT solutions generically possess conical 
singularities. In order to avoid such singularities, the following function must go to unity 
on the symmetry axis 
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where X is the squared norm of the axisymmetric Kilhng vector. In Weyl coordinates 
this function behaves hke e~^°" as p — 0. The condition (28) is thus equivalent to the 
requirement that the function a tend to zero in this hmit. Since a is invariant under the 
duahty transformation (11), one can simply compute it for the multi-black hole solutions. 
One gets, as a direct generalisation of the case of two positive mass black holes given in 
[12], that 

1=1 1<J ^ J ^ J 



where 



Ef^ = ri±rj± + {z - Zi± Ci){z - Zj ± Cj) + . (30) 



On the segment Dk, the function a is thus constant and is equal to 



k-l h 



\ ^ \ ^ I {Zi — Zj)"^ — {Ci -\- c-^"^ 



i=l j=k ^J')^ 

= sign (qc.) In I 1 - ^ ) (31) 

where La = \zj Zj I I C-j I I Cj I IS the distance between the two horizons of the black 
holes centred aX z = Zi and z = zj respectively. Since we require the horizons not to 
overlap, all the Lij are strictly positive and one sees that a can only be zero on each 
segment if some of the masses Cj are negative. 

Our multi-NUT solution defines thus a perfectly smooth space-time outside the hori- 
zons if and only if 

Yitl f'~%~f'^% = l (32) 

for all k between 2 and h. These h — 1 equations determine the relative positions of the 
NUT sources as functions of their charges. 



4 Some examples 

Let us consider a simple class of examples with three NUT sources, with charges rii = 
Us = puQ and n2 = —qriQ for two integers p and q. We also fix z^ = —Z2 = Zq and Zi = 0. 
The absence of a conical singularity requires that 

(^o^-(p-g)W)((2^o)'-(2pno)2) _^ 



(V-(p + g)W)(2^o) 



(33) 
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This equation can be solved for p > 4g, by 



p-q 



1 _ i£ 

P 



-.no 



(34) 



and the horizons are disjoint for any value of p and q. 

The asymptotic r = const slices are then diffeomorphic to a Lens space S^fE^ with 
N = 2p — q. The value of the Chern class in these examples is = 2r + 7g for strictly 
positive integers r and q. 

Since ^Bdip defines the pullback of a U{1) connection on each Ui of V, we can com- 
pute the Chern class of the two-cycles in V from it. We define the partitions of various 
two-cycles over the atlas of V as depicted in the following figures 




Partition on the open set Ui 




Partition on the open set C/3. 




Partition on the open set 1/2. 




Partition on the open set C/4. 



The Chern class of the cycle S*" is given by 

No. = / dB\u,Adif+-^ / dB\u,Adip = / {B\u, - B\u,)d^ = p (35) 

87^0 Js§ STmo Js-^ 87ino Jgg^ 

and one computes in the same way that A''^ = —q, = p — q and Ng = 2p — q. 
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5 Gravitational instantons 



We recall that the quantum mechanics of a particle in a Taub-NUT space-time requires 
the quantisation of the product of its mass with the NUT charge of space-time [13], ex- 
actly as in the case of a magnetic monopole. Moreover, the NUT charge as defined in [14] 
for not necessarily stationary space-times is shown to be preserved by small deformations 
of the solutions through the introduction of gravitational waves. 

The Euclidean self-dual Taub-NUT solutions might play a role in quantum gravity 
very similar to the one played by instantons in gauge theories [15]. The analogue of the 
instanton number would then be given by the Chern class of the asymptotic Lens space, 
in the sense that the action evaluated for such a solution is proportional to |A^|. The 
index of the Dirac operator is however given by the Pontryagin number. 

The solitons we have described in this letter are the Minkowski analogues of the 
instantons described in [5] with the slight generalisation of considering both positive 
and negative mass. However, the singularities associated with negative masses are not 
removed by the effects of the NUT charges in the Euclidean case. 

One can Wick rotate the Minkowskian solitons to Euclidean-signature solutions by 
choosing a complex pure imaginary parameter for the duality transformation (20). In- 
deed, for the Riemannian metric in Weyl coordinates, 

ds^ = H-^e^" [dz^ + dp^) + p^H^^d^^ + H{d^ + Bd^pf (36) 

the Euclidean Ernst equation is 




where the real Ernst potentials are S± = H ± B, with B derived from B using equation 
(15). For a static Ernst potential, i.e. one satisfying £"+ = the Euclidean Ernst 
equation is identical to the Minkowski one, and the multi-black hole solutions are thus 
solutions of the Euclidean theory as well. The Euclidean Ernst equation is the equation 
of motion of an SL{2, R) / 5*0(1, 1) non-linear sigma model, and it is left invariant by the 
5*0(1, 1) Ehlers transformation 

cosh a £:± =F sinh a 

£±{a) = r r-r — ^ • (38) 

cosh a =F smh a t± 

Applying this transformation, we obtain the following potentials for the Riemannian 
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metric (36)^ 



H = ^— \ ^ B = h-2yn,^-^ (39) 

cosh a \\{ri + q)^ — sinh a ~ q)^ ~^ 

where is defined as in the Minkowski case and the mass and the NUT charges are 
given by mj = cosh 2a q and rij = sinh 2a Cj. The resolution of the Dirac-Misner string 
singularities goes the same way. All the NUT charges are thus required to be integral 
multiples of a fundamental NUT charge no, and the imaginary time coordinate ip is again 
periodic, with period Svrno. 

However one can not get rid of the conical singularities in the Euclidean case without 
introducing singularities associated with negative masses. The only regular instantons 
left over are thus the single instanton with m = ||n| and the self-dual instantons for 
which Cj = [5]. 

The (anti) self-dual gravitational instantons with = irij [16] can be obtained by 
taking the limit Cj — > 0, a ^ ±oo while holding cosh 2a Cj fixed and equal to mj. In this 
limit, the Ernst potentials behave as 

I — g±2a Ci_ 

^± = i + + ^(^^') (40) 



and one computes that the function rj becomes ^y p'^ + {z — Zi)'^ and 



1=1 



The Ernst potentials then verify Szp = 1 for n.i = ±mi respectively, and the Ernst equation 
reduces to the linear differential equation 



(42) 



For (anti) self-dual instantons Szp = const, cr = and equation (32) turns out to be 
satisfied independently of the position of the sources on the axis. However, the absence 
of Euclidean NUT singularities nevertheless requires all masses to be equal to no. 
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